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Abstract 



We have made a precise measurement of the central inclusive jet cross section 
at y/s = 1.8 TeV. The measurement is based on an integrated luminosity of 
92 pb^ 1 collected at the Fermilab Tevatron pp Collider with the D0 detector. 
The cross section, reported as a function of jet transverse energy (Ej- > 
60 GeV) in the pseudorapidity interval \r/\ < 0.5, is in good agreement with 
predictions from next-to-leading order quantum chromodynamics. 
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Typeset using REVTpjX 



Within the framework of quantum chromodynamics (QCD), inelastic scattering between 
a proton and an antiproton can be described as an elastic collision between a single proton 
constituent and a single antiproton constituent. These constituents are often referred to as 
partons. After the collision, the outgoing partons manifest themselves as localized streams 
of particles or "jets". Predictions for the inclusive jet cross section are given by the folding 
of parton scattering cross sections with experimentally determined parton distribution func- 
tions (pdf's). These predictions have recently improved with next-to-leading order (NLO) 
QCD scattering calculations [1-3] and new, accurately measured pdf's [4,5]. We measure the 
cross section for the production of jets as a function of the jet energy in the plane transverse 
to the incident beams, Et- The measurement is based on an integrated luminosity of 92 
pb^ 1 of pp collisions collected with the D0 detector || at the Fermilab Tevatron Collider. 
Measurements of inclusive jet production with smaller integrated luminosity have been per- 
formed previously by the UA2 and CDF collaborations [7,8]. The cross section measurement 
presented here allows a stringent test of QCD, with a total uncertainty substantially reduced 
relative to previous results. 

Jet detection in the D0 detector utilizes primarily the uranium-liquid argon calorimeters 
which have full coverage for pseudorapidity \rj\ < 4.1 (rj = — ln(tan(#/2)), where 9 is the 
polar angle relative to the proton beam). The calorimeter is segmented into towers of 
A?7 x A<f) = 0.1 x 0.1, where <ft is the azimuthal angle. 

Initial event selection occurred in two hardware trigger stages and a software stage. 
The first hardware trigger selected an inelastic pp collision - indicated by signals from the 
trigger hodoscopes located near the beams on either side of the interaction region. The 
next stage required transverse energy above a preset threshold in calorimeter trigger tiles 
of At] x A0 = 0.8 x 1.6. Selected events were digitized and sent to an array of proces- 
sors. Jet candidates were then reconstructed with a cone algorithm and the entire event 
recorded if any jet Et exceeded a specified threshold. For software jet thresholds of 30, 50, 
85, and 115 GeV, integrated luminosities of 0.34, 4.6, 55, and 92 pb^ 1 , respectively, were 
accumulated during a 1994-1995 data run. The two lowest Et triggers required a single 
proton-antiproton interaction within the beam crossing as signaled by timing information 
in the trigger hodoscopes. 

Jets were reconstructed offline using an iterative fixed-cone algorithm with a cone radius 
of 71=0.7 in space [[|. Background from isolated noisy calorimeter cells and accelerator 
beam losses which mimicked jets were eliminated with quality cuts |L(J. Background events 
from cosmic ray bremsstrahlung or misvertexed events were eliminated by requiring the 
missing transverse energy in each event to be less than the larger of 30 GeV or 0.3i?™ ax , 
where i?™ ax is the Et of the leading jet. Residual jet contamination is less than 1% at all 
Et, based on event simulations with superimposed calorimeter noise distributions and on 
visual scanning of jet candidates with Et greater than 350 GeV. The jet selection efficiency 
for \t]\ < 0.7 has been measured as a function of jet E T and found to be (97±1)% below 
250 GeV and (95±2)% at 400 GeV. 

At high instantaneous luminosity, more than one interaction in a single beam crossing is 
probable ( ~ 20% for this data set). The event vertex was reconstructed using data from the 
central tracking system. For events with multiple vertices, the two vertices with the largest 
number of tracks were retained. Due to fluctuations of jet charged-particle multiplicity, an 
additional parameter was used to select the vertex. If an event had more than one vertex, 
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the quantity St = \EE^\ was calculated for both vertices. The vertex with the smaller St 
was selected as the event vertex and used to calculate jet Et and 77. The selected vertex was 
required to be within 50 cm of the detector center. This last requirement retained (90±1)% 
of the events, independent of jet Et- 

The transverse energy of each jet was corrected for the underlying event, additional 
interactions, noise from uranium decay, the fraction of particle energy showered outside of 
the reconstruction cone, detector uniformity, and detector hadronic response. A complete 
discussion of the jet energy scale calibration can be found in Ref. ||11|| . For \r]\ < 0.5, the mean 
total correction factor for jet energy is 1.154±0.017 [1.118±0.023] at 100 GeV [400 GeV]. 

The inclusive jet cross section was computed in contiguous ranges of Et using data from 
the four trigger sets. The spectrum includes data from the 30 GeV trigger between 60- 
90 GeV, from the 50 GeV Et trigger between 90-130 GeV, then from the 85 GeV trigger 
between 130-170 GeV, and above 170 GeV from the 115 GeV trigger. The single interaction 
requirement on the two lowest-.&r triggers introduced an inefficiency which was corrected 
by matching the 50 GeV trigger cross section to the 85 GeV trigger cross section above 130 
GeV, where both triggers are fully efficient. This introduces an additional 1.1% luminosity 
uncertainty to the 50 GeV trigger set. A similar matching between the lowest-Er trigger and 
the 50 GeV trigger introduces another 1.4% uncertainty for the lower set, which is added in 
quadrature to the 1.1% matching uncertainty. 

The steep E T spectrum is distorted by jet energy resolution. At all E T , the resolution 
(measured by balancing Et in jet events) is well described by a gaussian distribution; at 
100 GeV the standard deviation is 7 GeV. The distortion was corrected by assuming an 
ansatz function (AEt~ b ) ■ (1 — 2Et / \/s) c , smearing it with the measured resolution, and 
comparing the smeared result with the measured cross section. The procedure was repeated 
by varying parameters A, B, and C until the best fit was found between the observed cross 
section and the smeared trial spectrum. The ratio of the initial ansatz to the smeared ansatz 



was used to correct the cross section on a bin- by-bin basis ||12|| . The resolution correction 
reduces the observed cross section by (13±3)% [(8±2)%] at 60 GeV [400 GeV]. 

The resulting inclusive jet cross section for < 0.5, shown in Fig. [I], has been averaged 
over each Et bin (AEt) and over the central unit of rapidity (A?7=l). This bin-averaged 
double differential cross section, (d 2 a/(dE T dr])), was calculated as N/( AE T ArjeC) where N 
is the total number of jets observed in the bin, e is the selection efficiency as a function of 
Et, and £ the integrated luminosity associated with the trigger set. The cross section is 



consistent with a preliminary measurement from a smaller 1992-1993 data set [10| . 

Figure [I] also shows a theoretical prediction for the cross section from the NLO event 
generator JETRAD ||. There is good agreement over seven orders of magnitude. Inputs to 
the NLO calculation are the renormalization scale \x (always chosen to equal the factorization 
scale), the parton distribution function (pdf), and the parton clustering algorithm. For the 
calculation shown here, \l = 0.5-E'™ ax and the pdf is CTEQ3M [[|. Partons separated by 
less than lZ mp =1.31Z were clustered if they were also within 71=0.7 of their .^-weighted 
r]-(p centroid. This choice of 7^ sep is discussed in Ref. 0. Variations in the predicted cross 
section due to the input choices are approximately 30% fi~3|j . 

The data in Fig. [I] have an overall luminosity uncertainty of 6.1%. The data are plotted at 
the Et value for which a smooth function describing the cross section is equal to the average 
cross section in the bin. The band shows the total systematic uncertainty as a function of 
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The \rj\ < 0.5 cross section. 




Plotted E T 


Bin Range 


Cross Sec. ± Stat. 


Sys. 


(GeV) 


(GeV) 


(fb/GeV) 


Uncer.(%) 


64.6 


60- 70 


(6.59±0.04)xl0 6 


±8 


74.6 


70- 80 


(2.89±0.03)xl0 6 


±8 


84.7 


80- 90 


(1.41±0.02)xl0 6 


±8 


94.7 


90 - 100 


(7.07±0.04)xl0 5 


±8 


104.7 


100 - 110 


(3.88±0.03)xl0 5 


±8 


114.8 


110 - 120 


(2.21±0.02)xl0 5 


±8 


124.8 


120 - 130 


(1.27±0.02)xl0 5 


±8 


134.8 


130 - 140 


(7.70±0.04)xl0 4 


±8 


144.8 


140 - 150 


(4.86±0.03)xl0 4 


±8 


154.8 


150 - 160 


(3.07±0.02)xl0 4 


+9, - 8 


164.8 


160 - 170 


(2.00±0.02)xl0 4 


±9 


174.8 


170 - 180 


(1.34±0.01)xl0 4 


±9 


184.8 


180 - 190 


(9.12±0.10)xl0 3 


±9 


194.8 


190 - 200 


(6.15±0.09)xl0 3 


+10, - 9 


204.8 


200 - 210 


(4.29±0.07)xl0 3 


±10 


214.8 


210 - 220 


(2.93±0.06)xl0 3 


+11, -10 


224.8 


220 - 230 


(2.14±0.05)xl0 3 


+11, -10 


239.4 


230 - 250 


(1.30±0.03)xl0 3 


±11 


259.4 


250 - 270 


(6.54±0.20)xl0 2 


+12, -11 


279.5 


270 - 290 


(3.77±0.15)xl0 2 


+13, -12 


303.9 


290 - 320 


(1.79±0.08)xl0 2 


+15, -13 


333.9 


320 - 350 


(6.82±0.52)xl0 1 


+17, -15 


375.7 


350 - 410 


(1.89±0.19)xl0 1 


+20, -17 


461.1 


410 - 560 


(1.24±0.31)xl0° 


+30, -26 



E T . Listed in Table | are the plotted values of E T , the E T ranges, the cross section, and the 
statistical and systematic uncertainty. The tabulated systematic uncertainties include jet 
and event selection, unsmearing, relative luminosity, and energy scale uncertainties added 
in quadrature. The 6.1% luminosity uncertainty is not included. 

Figure H shows the various uncertainties for the \rj\ < 0.5 cross section. Each curve 
represents the average of the nearly symmetric upper and lower uncertainties. The uncer- 
tainty attributed to the energy scale varies from 8% at low E T to 30% at 450 GeV. This 
contribution dominates all other sources of uncertainty, except at low E?, where the 6.1% 
luminosity uncertainty is of comparable magnitude. 

The |?7| < 0.5 region provides our optimum test for departures of data from NLO QCD. 
In this region, the detector is uniformly thick (seven or more interaction lengths with no 
gaps) and both jet resolution and calibration are precise. Also, jet production from the 
scattering of possible constituents within quarks is largest for 77 = 0, relative to standard 
QCD predictions (T^|. For comparison to Ref. ||, we have also carried out a similar analysis 
in the region 0.1 < \r]\ < 0.7. 

Figure | shows the ratios (D - T) /T for the data (D) and JETRAD NLO theoretical (T) 
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TABLE II. Correlations of the total uncertainty in the cross section. 



£r(GeV) 64.6 104.7 204.8 303.9 461.1 

646 LOO 096 085 071 O40~ 

104.7 096 1.00 092 0.79 0.46 

204.8 0.85 092 1.00 0.91 0.61 

303.9 0.71 0.79 0.91 1.00 0.67 
461.1 0.40 0.46 0.61 0.67 1.00 



predictions based on the CTEQ3M, CTEQ4M and MRST pdf's [4,5] for \rj\ < 0.5. Given 
the experimental and theoretical uncertainties, the prediction is in good agreement with 
the data; in particular, the data above 350 GeV show no indication of an excess relative to 
QCD. 

The data and theory can be compared quantitatively with a x 2 = Y^i,j{Di — 
Ti)(C~ 1 )ij(Dj — Tj) test incorporating the uncertainty covariance matrix C . Here Di and 
Ti represent the z-th data and theory points, respectively. The elements CV,- are constructed 
by analyzing the mutual correlation of the uncertainties in Fig. ||] at each pair of Et values. 
As shown, the luminosity and resolution uncertainties are quite significant at low Et and 
uncertainties in the jet energy scale dominate at high Et- Since the luminosity uncertainties 
are fully correlated over all Et, as are the resolution uncertainties, and since the energy scale 
uncertainties at large Et are also correlated, the overall systematic uncertainty is largely 
correlated. Table [TI] shows that the bin-to-bin correlations in the full uncertainty for rep- 
resentative Et bins are greater than 40% and positive. (The full covariance matrix can be 
found in Ref. 0.) 

Table |1| lists x 2 values for several JETRAD predictions incorporating various parton 
distribution functions [4,5]. Each comparison has 24 degrees of freedom. The JETRAD 
predictions have been fit to a smooth function of Et- All five predictions describe the 
\rj\ < 0.5 cross section very well (the probabilities for x 2 to exceed the listed values are 
between 47 and 90%). The 0.1 < \r]\ < 0.7 cross section is also well described (probabilities 
between 24 and 72%). Comparisons for \rj\ < 0.5 with EKS [ffl| calculations using CTEQ3M, 
TZ sep =1.3TZ, and with renormalization scales \x = cE™ 3 "* or cE 1 ^ 1 , where c varies, are shown in 
Table [TV]. These calculations also describe the data very well (better than 57% probability) 
at all renormalization scales. 

The top panel in Fig. ^ shows (D — T)/T for our data in the 0.1 < \rj\ < 0.7 region 
relative to an EKS calculation using the CTEQ3M pdf, /i = 0.5E^\ and n sep =2.01l. (The 
tabulated data can be found in Ref. [|Lj.) Also shown are the data of Ref. pi relative to 
the same EKS prediction. For this rapidity region, we have carried out a x 2 comparison 
between our data and the nominal curve describing the central values of the data of Ref. || . 
Comparing our data to the nominal curve, as though it were theory, we obtain a x 2 of 63.2 
for 24 degrees of freedom (probability of 0.002%). Thus our data cannot be described with 
this parameterization. As illustrated in the bottom panel of Fig. |], our data and the curve 
differ at low and high Et] such differences cannot be accommodated by the highly correlated 
uncertainties of our data. If we include the systematic uncertainties of the data of Ref. || 
in the covariance matrix, the x 2 is reduced to 24.7 (probability of 46%). 

In conclusion, we have made the most precise measurement to date of the inclusive jet 
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TABLE III. x 2 comparisons between JETRAD and \rj\ < 0.5 and 0.1 < \rj\ < 0.7 data for 
fj, = Q.bEtp 8 *, 7£ sep =1.37*!., and various pdfs. There are 24 degrees of freedom. 



pdf \rj\ < 0.5 0.1 < jgj < 0.7 

CTEQ3M 23^9 28^4 

CTEQ4M 17.6 23.3 

CTEQ4HJ 15.7 20.5 

MRSA' 20.0 27.8 

MRST 17.0 19.5 



TABLE IV. x 2 comparisons between EKS and \rj\ < 0.5 data for CTEQ3M, ft sep =1.37£, and 
various renormalization scales given by cEf ax or cE^ 1 . There are 24 degrees of freedom. 



c 


...2 /'iTimaxN 


x 2 (4 et ) 


0.25 


14.8 


19.8 


0.50 


19.4 


22.2 


1.00 


16.8 


18.1 



cross section for Et > 60 GeV. QCD predictions are in good agreement with the observed 
cross section for standard parton distribution functions and different renormalization scales. 



This is consistent with our previous measurements of dijet angular distributions [ 141 , which 



are also in good agreement with QCD and show no indication of physics beyond the Standard 
Model. 
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FIG. 1. The \r/\ < 0.5 inclusive cross section. Statistical uncertainties are invisible on this scale. 
The solid curves represent the ±la systematic uncertainty band. 
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FIG. 2. Contributions to the |r/| < 0.5 cross section uncertainty plotted by component. 
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FIG. 3. The difference between data and JETRAD QCD predictions normalized to predictions. 
The bands represent the total experimental uncertainty. 
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FIG. 4. Top: Normalized comparisons of our data to EKS and of the data in Ref. [|| to EKS. 
Bottom: Difference between our data and smoothed results of Ref. [Q normalized to the latter. 
The band represents the uncertainty on our data. 
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